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Abstract. In this paper, we present a microscopic model for heterogeneous ferroelectric and an 
order parameter for relaxor phase. We write a Landau theory based on this model and its application 
to ferroelectric PbFei/2Tai/203 (PFT) and antiferroelectric NaNbOsiGd. We later discuss the 
coupling between soft mode and domain walls, soft mode and quasi-local vibration and resulting 
susceptibility function. 



A simple model ferroelectric thin film with dead layers was recently analyzed by 
Bratkovsky and Levanyuk [1] resulting in an analytic solution. Based on this study, we 
proposed a microscopic model for an inhomogeneous ferroelectric consisting of ferro- 
electric slabs sandwiched between dielectric interfacial layers [2]. Low energy solutions 
of these models reveal that the ferroelectric slabs break into alternating domains (Fig. 1) 
with zero total macroscopic polarization. The size of the domains was shown to depend 
only on the relative total width of the dielectric and ferroelectric regions in the direction 
of the field [2]. The nanodomain structure appears cooperatively and its origin lies in the 
reduction of depolarization field [1,2]. 

The alternating polarization domains are accompanied by shear strain owing to the 
electro strictive coupling. For the case of domains alternating along (110) direction, we 
obtained atomic displacements corresponding to polarization and strain fields and calcu- 
lated diffused scattering intensity [2], which agrees well with recent neutron scattering 
data [3]. We suggest that the relaxor phase corresponds to a finite value of the order 
parameter describing these alternating domains. 

The fluctuations of the order parameter introduced above are conjugated with the field 
which we will call Hq. We will assume the existence of quenched field Hq below Burns 
temperature. We use the notation rj for the relaxor order parameter (the magnitude of 
polarization inside a nanodomain) and write the Landau expansion: 



F = Fo + l-aP^ + IpP^ + \yP^ -f 
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where a = ao(r — Tciy) and A = a(r — T^). 

At zero macroscopic polarization, the equilibrium condition with respect to r\ is 
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FIGURE 1. A model of inhomogeneous ferroelectrics [2]. 
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At high temperatures r\ vanishes as (T — 7,^)^^ In this limit, dielectric permittivity 
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there is a deviation from this law, and the dielectric 



obeys a Curie-Weiss law. At T 
permittivity peak is diffused due to the coupling between the new order parameter 
appearing at the phase transition and frozen conjugated fields. 

We used the expression derived in order to describe the diffuseness of the phase 
transition in relaxor PFT from the para-phase to the relaxor phase with ferroelectric 
nano-regions. We also took into account a low temperature (glass-type) phase transition 
which results in a strong decrease of dielectric permittivity below approximately 220 K. 
We fitted expressions obtained to experimental data (Fig. 2) and got the best fit shown 
by the solid line. We find that is rather close to the extrapolated high-temperature 
Curie-Weiss temperature. This justifies that the ferroelectric relaxor order parameter, 
T], is connected with local polarization, and can be regarded to a wave of polarization 
described above. 

We use the same expression (3) to treat the diffused phase transition in an antiferro- 
electric NaNbOs doped with Gd, where the significance of r\ now is an antiferrolectric 
order parameter. The resulting fit to experimental data is rather good shown in Fig. 3. 
The critical temperatures are indicated in the figure. We also show in Fig. 3 the differ- 
ence between the inverse dielectric permittivity and high temperature Curie-Weiss law. 




It is seen that this difference, at low temperatures, behaves linearly with temperature 
while at the dielectric permittivity peak position it is diffused in accord with the theory. 

Experimental data [4] show that the "waterfall phenomenon" exists in the whole 
temperature interval between the freezing temperature and Bums temperature, implying 
that the soft mode is strongly damped in the relaxor phase (note that this soft mode 
is not the uniform ferroelectric). We believe this mechanism to be connected with the 
interaction of the soft mode vibrations with domain walls and with local dipoles (a 
mathematical expression is similar to that obtained for the soft mode coupled with 
microscopic dipoles [5, 6]). 

The acoustic mode can have a dip due to its coupling with optical mode [7] (Fig. 4a). 
We suggest that, due to disorder, the deviation of the optical and acoustic modes from 
average is different in different regions of the disordered crystal and one has to introduce 
a distribution function of these deviations (Fig. 4b). We also consider the case when the 
optical mode has a dip (Fig. 4c) if the interaction with the acoustic mode is not taken 
into account (the first row in Fig. 4). In this case the optical mode dispersion curve is: 
e = (Dq — kq^ + bq^ + .... We derived in this case that the correlation function deviates 
from the Omstein-Cemike expression by an oscillating factor: 

X{r)^^ sin(/ro^) exp(-/ccr) (4) 

where kc is inverse correlation length and is the wave vector of the susceptibility 
oscillations. In the ^-space susceptibility looks as 




FIGURE 3. Temperature dependencies of e' (1) and lO'^/e' (2) measured at 10^ Hz for O.SSNaNbOs- 
0.12Gdi/3NbO3 crystal, and the difference between the experimental 1/e' dependence and the Curie- 
Weiss fit of 1/e' (3). Solid lines are the best fits of the theory to experimental data. 
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Here k^ = kQ — and — Ikokc- These constants are expressed in terms of parameters 
in the optical mode dispersion: (Oq /b = k^ + and k/b = 2k^. The inverse correlation 
length kc can be obtained from: 
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where, at Oq > ^ > 0, kc is real. Otherwise, it is imaginary, and there appear harmonic 
beatings. 

A general form of the Hamiltonian taking into account the interaction between the 
acoustic (m) and optical (x) displacements (and a dip in the optical mode) is: 



Hharmiq) = ^U-qA{q)Uq + U-qVXq +X-qV*Uq + ^X-qB{k)Xq 
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where, at small q. 



A{q) = aq^ 



(8) 
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FIGURE 4. Model dispersion curves: a) the condensation of the modes is due to mode-mode coupling, 
b) the same as (a) but with taking into account disorder, c) the condensation of the modes is due to both 
the polarization instability and mode-mode coupling; The first row is without the mode-mode coupling, 
the second row is with the mode-mode coupling; The dashed region arises because of disorder. 

B{k) = + (9) 

V{q) = Viq^ + iv2Pq (10) 

Here P is polarization; The imaginary part of the coupling term appears in symmetry 
broken regions and in the regions where there exists gradient of polarization, that is at 
the boundaries of the domains or/and polar regions. This coupling leads to a "repulsion" 
of the acoustic and optical modes as it is shown in Fig. 4. 

We also will discuss the correlation of polarization in the case of a finite volume of 
the polar region. Let us write dielectric permittivity of a finite uniform volume V: 

X-^J dVdV'{P{r)P{r')) ~ ^mm(r2,r2) (11) 

where an Omstein-Cernike correlator was used, K being the constant in front of 
(VP) in the Landau fluctuation energy, and rc ~ (r — Tew) > iii the first approx- 
imation. It is seen from this expression that, at high temperatures, susceptibility behaves 
according to the Curie- Weiss law; when the ferroelectric correlation radius reaches the 
volume V size then susceptibility saturates. The decrease in polarization at the boundary 
of a polar region would smoothen this crossover. 

We finally consider the Hamiltonian consisting of the part describing the long-range 
ordered polarization Pq (with a wave vector q) and coupled with it local polar vibrations 
with local polarization P/: 



+ \bP^ + ^gP^Pi - EP^ - EPi (12) 



The susceptibility connected with quasilocal vibrations at zero polarization Pq can be 
found by an ordinary procedure [5]: 



xi = YZ 

where 
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It is seen that the quasilocal vibrations become unstable when OCqiXqO < 1 • It implies 
that the local vibrations satisfying this condition will freeze in and will be arranged in 
space in accordance with the wave vector q. These results are consistent with resent 
experimental findings [8] that there is a local transformation at Bums temperature. 
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